library(tidyverse)
data_set <- read_rds("../data/aula6-01.rds")
glimpse(data_set)
## Rows: 45
## Columns: 10
## $ macrofita <chr> "BRASU", "BRASU", "BRASU", "BRASU", "BRASU", "BRASU", "BRASUâŠ
## $ dose_t_ha <dbl> 13.5, 13.5, 13.5, 13.5, 13.5, 27.0, 27.0, 27.0, 27.0, 27.0, âŠ
## $ repeticao <dbl> 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, 1, 2, 3, 4, 5, âŠ
## $ mo <dbl> 8, 7, 8, 7, 7, 10, 10, 9, 9, 9, 9, 10, 10, 8, 10, 8, 8, 8, 7âŠ
## $ p_h <dbl> 5.8, 5.8, 5.9, 6.0, 6.0, 6.1, 6.0, 5.7, 6.1, 6.1, 6.1, 6.1, âŠ
## $ p <dbl> 3, 4, 4, 3, 4, 4, 4, 3, 5, 4, 4, 3, 3, 3, 3, 8, 7, 7, 7, 10,âŠ
## $ k <dbl> 1.01, 0.86, 5.64, 4.57, 4.71, 11.86, 11.95, 7.78, 9.33, 7.70âŠ
## $ ca <dbl> 17.95, 19.95, 18.98, 15.81, 16.37, 21.79, 21.03, 17.19, 20.6âŠ
## $ mg <dbl> 6.29, 6.90, 6.64, 5.57, 5.76, 7.50, 7.55, 6.14, 6.91, 5.93, âŠ
## $ al <dbl> 0.28, 0.25, 0.19, 0.26, 0.20, 0.21, 0.22, 0.22, 0.20, 0.14, âŠ
library(DT)
datatable(data_set)
DICA:PARA COLOCAR NEGRITO, COLOQUE A PALAVRA ENTRE 4 ASTERISCOS
DICA: PARA COLOCAR EM ITĂLICO, COLOQUE A PALAVRA ENTRE DOIS ASTERISCOS
variaveis <- data_set %>%
select(mo:al) %>%
names()
variaveis
## [1] "mo" "p_h" "p" "k" "ca" "mg" "al"
for(i in seq_along(variaveis)){
my_var <- variaveis[i]
macrofitas <- data_set %>% pull(macrofita)
doses <- data_set %>% pull(dose_t_ha)
y <- data_set %>% pull(my_var)
print("------------------")
print(my_var)
print("------------------")
media <- mean(y,na.rm=TRUE)
desv_pad <- sd(y,na.rm = TRUE)
cat("Média:",round(media,2),"\n")
cat("Desvio-padrĂŁo:",round(desv_pad,2),"\n")
# GrĂĄfico de histograma
plot_hist <- tibble(macrofitas,doses,y) %>%
ggplot(aes(x=y))+
geom_histogram(bins = 7,fill="grey",colour="black")+
theme_bw()+
labs(x=my_var)
print(plot_hist)
# AnĂĄlise de variĂąncia
ExpDes.pt::fat2.dic(macrofitas,doses,y,
sigF = 0.15,sigT = 0.05,
mcomp = "tukey",
quali = c(TRUE,FALSE),
fac.names = c("macro","dose"))
}
\[ y_{ij}=\mu + \tau_{i} + \beta_{j} + (\tau \beta_{ij})_{k} + \epsilon_{ijk} \]
| FV | GL |
|---|---|
| A | \((a-1)\) |
| B | \((b-1)\) |
| AB | \((a-1)(b-1)\) |
| R | \(ab(r-1)\) |
| T | \(abr-1\) |
library(plotly)
graf <- data_set %>%
group_by(dose_t_ha,macrofita) %>%
summarise(mo=mean(mo)) %>%
ggplot(aes(x=dose_t_ha,y=mo,fill = macrofita))+
geom_col(position = "dodge",colour="black")+
theme_bw()+
scale_fill_viridis_d()
ggplotly(graf)
## Typing cat (meme)